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INTRODUCTION 

A  rigorous  ray  tracing  theory  has  already  been  developed  to  first  order 
where  it  was  shown  that  an  effective  sound  velocity  profile  was  in  general 

dependent  upon  frequency.   A  higher  order  theory  for  wave  mechanics  has 

2 
also  been  developed  in  one  dimension.   Exact  albeit  numerical  wave  functions 

have  recently  been  developed  for  bound  states  (i.e.,  trapped  modes  in  under- 
water acoustics)  in  one  dimension  where  the  turning  point  (vertex  depth  in 
underwater  acoustics)  of  the  phase  integral  have  been  shown  to  recede  to  °° 
for  indices  of  refraction  that  remain  finite  and  where  the  effective  action 
variable  (the  phase  integral  around  the  branch  cut  between  the  turning  points 
for  the  N   eigenvalue,  N  =  1,2,3,  ...)  in  units  of  phase  has  been  shown  to 

be  N  cycles  for  the  N   trapped  mode.   The  modes  have  been  shown  to  be 

3  4 
numerically  exact  within  the  limits  of  computer  round-off.  * 

In  this  exposition,  the  terminology  "rigorous  ray  tracing"  denotes  an 
exact,  albeit  numerical  ray  tracing  while  "classical  ray  tracing"  denotes 
the  conventional  ray  tracing  that  is  valid  in  the  short  wavelength  or 
infinite  frequency  limit. 

Herein  we  apply  some  of  these  developments  in  wave  mechanics  to  under- 
water acoustics  in  order  to  generate  for  depth-dependent  sound  velocity  pro- 
files a  rigorous  ray  theory  that  is  exact,  albeit  numerical,  and  that  is 
mutually  consistent  with  a  normal  mode  theory.  The  objectives  of  this  inves- 
tigation are  to  gain  insight  into  the  physics  of  rigorous  ray  tracing  from 
a  few  examples  of  the  exact  effective  sound  velocity  profile,  to  develop  the 
correct  equations  of  motion  for  the  ray,  and  to  show  that  the  ray  directions  in 
general  are  not  orthonormal  to  the  wave  surfaces  of  constant .phase  (i.e., 


the  wave  normals  and  the  ray  directions  for  the  same  points  are   not  in 
general  paral lei ) . 

In  this  exposition  we  assume  that  the  sound  velocity  profile  is  a  func- 
tion of  depth  only  and  that  for  finite  depths,  the  inverse  square  of  the 
sound  velocity  profile  is  finite.  The  surface  is  at  zero  depth.  The  earth 
is  assumed  to  be  flat  and  we  assume  a  bottomless  ocean.  Herein  cylindrical 
coordinates  are  used  and  we  assume  the  source  is  located  on  the  z-axis  at 
some  depth,  z  .  Hence  we  have  azimuthal  symmetry.  We  also  assume  that  the 
sound  source  is  a  point  harmonic  source  whose  frequency  is  oo  in  radians  per 
second.  All  numerical  computations  were  based  upon  a  Runge-Kutta  fourth- 
order  technique  with  a  step  size  of  25  m. 

In  Section  I  of  this  exposition,  we  develop  Hamilton's  principal  func- 
tion, S,  for  rigorous  ray  tracing  which  is  based  upon  the  alternative  de- 
velopment presented  in  Appendix  A  of  Ref.  1.  (Heretofore  S  had  been  cited 
as  the  effective  Hamilton's  principal  function.)  Subsequently  we  develop 
the  equations  of  motion  by  the  Hamil ton-Jacobi  process  and  we  develop  the 
relationship  between  the  ray  direction  and  wave  normal.  In  Section  II  we 
generate  some  sample  exact  effective  sound  velocity  profiles,  C  (z)  from 
a  simple  sound  velocity  profile,  C(z).  From  these  exact  samples,  we  glean 
some  insight  into  the  behavior  and  application  of  rigorous  ray  tracing.  A 
brief  outline  of  applicable  portions  of  classical  ray  tracing  is  presented 
in  Appendix  A. 


I.   THEORY 

Here  we  shall  develop  the  equations  of  motions  for  a  ray   and  show  that 
the  ray  is  in  general  not  normal  to  the  acoustical  wave  surfaces  of  constant 
phase.  Here  for  convenience  and  without  loss  of  generality,  we  express  the 
sound  velocity  profile,  C,  which  has  been  assumed  to  be  a  function  of  depth 
only,  by  the  form 


C(z)   =   [A  -  B(z)f1/2 


where  z  is  depth  and  where  A,  a  real  constant,  and  B(z),  a  real  function, 
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are  chosen  such  that  at  the  source  depth,  z  ,  of  the  sound,  C(z  )  =  A 

As  C(z)  is  only  depth  dependent,  the  Helmholtz  equation  is  separable 

such  that  the  z  component  may  be  expressed 


826(z)/8z2  +  uj2[A  -  C~2  -  B(z)]6(z)  =  0  (1) 


_2 
where  6  is  a  scalar,  oj  is  frequency  in  radians  per  second,  C   is  the  real 

_2 
separation  constant  such  that  A  >  C  .  In  addition,  C  is  also  the  vertex 

(WKB  turning  point)  velocity.  Let  us  use  an  ansatz,  which  has  been  inspired 

by  the  WKB  approximation,  for  6  of  the  form 


6±  =  (A  -  C2   -  e)"1/4  expfxi/  o,(A  -  C"2  -  B)1/2dz'].        (2) 


Substituting  Eq.    (2)   into  Eq.    (1)   generates  a  nonlinear  differential    equation 
for  our  new  variable,   £,   expressed  by 

6+     1      Iml +      5     /     3P/3Z      )2     .     B(z)  (3) 

4cu     (A-  Cc  -  p)       16c/ V  Cc  -  37 
m  m 


The  initial  conditions  for  0  require  that  Huygen's  construction  may  be 
accomplished  in  the  neighborhood  near  the  source.   Therefore  the  initial 
conditions  are 


■  B(z  ) 
s 
z  =  z 
s 


and 


93/3z|  =     dB/dz  | 

z=z  z=z 

s  s 

From  the   form  of  Eq.    (3),  we   see  that  6   is  a   function  of  z,  uj  and  C    .     We 

note  that   for  trapped  modes   Eq.    (3)   has  a  nodal-type  critical    pointr  such 

_2 
that  3  ->  A  -  C       from  below  as  z  ■*  °°.     This  node  at  z  =  °°  moves  the  vertex 
m 

depth  (turning  point)   for  the  ray   (for  flat  earth)    for  rigorous  ray  tracing 
from  the  WKB  vertex  depth  to   infinity.     Further   investigation  of  3  is  de- 
ferred to  Section  II . 

The  phase  of  our  ansatz   is  given  by 


/     co[A  -  3(z\u>,Cm)  -C~2]1/2dz' 


=     /     u)[C;2(z',u),Cm)-r2]1/2dz'     =     wW(z,u),Cm)  (4) 


where  by  analogy  to  Eq.    (A4)  U  is  Hamilton's  characteristic  function  for 
rigorous  ray  tracing  and  C     is  defined  to  be  the  effective  velocity  given 
by 


Ce(z,co,Cm)     =      [A  -   S(z,ui,Cm)f1/2.  (5) 


We  develop  Hamilton's   principal    function,    S,    for   rigorous    ray   tracing  as 


W(z,u,Cj    +  f  (6) 

m 


because  range  is  a  cyclic  coordinate  and  azimuth  is  symmetric.  By  Eqs .  (4) 

and  (6)  and  by  using  3S/3r  =  (1/C )  to  change  C  temporarily  from  a  function 

of  C  to  dS/dr,   we  have  the  Hamilton- Jacobi  equation  for  riqorous  ray 
m  ^  j 

tracing  expressed  by 

dS/dr  -    [C~e2(z,u,dS/dr)    -  (3S/9z)2]1/2  =  0.  (7) 

There  is  a  subtle  difference  between  Eqs.  (7)  and  (A2)  because  the  range 
dependence,  which  only  appears  explicitly  in  Eq.  (A2),  appears  in  Eq.  (7) 
both  explicitly  and  in  C  implicitly.  Thus,  even  though  C   is  a  separation 
constant  for  the  Helmholtz  equation,  C  is  no  longer  a  separation  constant 
for  Eq.  (7)  (albeit  Eq.  (6)  is  still  valid  for  cyclic  r) ,  but  (1/C  )  still 
is  a  transformed  momentum  of  the  Hamilton- Jacobi  process  for  both  Eqs.  (7) 
and  (A3).  Hence  C  is  in  general  anisotropic. 

Either  S  or  W  (at  the  worker's  option)  may  act  as  the  generating  func- 
tion from  which  the  equations  of  motion  for  rigorous  ray  tracing  may  be 
derived.  But  we  note  that  the  ansatz  may  be  expressed  as 


(8W/9z)"1/2  exp  ±iu)W. 


As  the  ansatz  describes  the  normal  mode  (the  analogous  exact,  albeit  numeri- 
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cal ,  quantum  mechanical  eigenfunctions  have  been  so  described  elsewhere  '  ), 

we  may  derive  the  eigenfunctions  from  W  by  quantizing  the  product  of  w  and 


the  action  variable  for  rigorous  ray   tracing  to  be  N  cycles  of  phase, 
N  =  1,2,3,  ...  .  Hence  W  is  the  keystone  by  which  we  may  derive  mutual ly 
consistent  normal  mode  and  rigorous  ray  tracing  theories  thus  rendering 
a  unifying  principle. 

The  wave  normals  for  rigorous  ray  tracing  are  given  by  the  gradients 
of  S  which  for  cylindrical  symmetry  renders 


VS  =  [C"2(z,ou,C  )  -C~2]1/2i  +  C_1i 
e      nr   m  J    z    m  r 


Thus  the  wave  normal  is  embedded  in  a  radial -vertical  plane  and  the  wave 
normal  may  be  described  by  the  deflection/elevation  (D/E)  angle,  \\>.     This 
D/E  angle  is  given  by 


C  C 

*  =  tan'1^)2-!]172  =  cos_1(^)  (8) 

e  m 


which  may  be  expressed  in  a  manner  similar  to  Snell's  law  for  classical 
ray  tracing  as 

C(z,u),Cj 

c   =  -s r^-  (8') 

m     cos  \\>{z) 

where  C  is  an  implicit  function  as  opposed  to  the  classical  ray  tracing, 
m 

Eq.  (Al),  where  it  is  an  explicit  function. 

We  now  derive  the  continuous  motion  for  rigorous  ray  tracing  for  the 
ray  passing  depth  z.  By  the  Hamil ton-Jacobi  process  the  range  r  that  a 
ray  passes  depth  z  is  given  for  cylindrical  symmetry  by 


as      n     aw 

3U7y  =      =  wy 


or 


z  C  [1  -  (C  /C  )3(&C  /3C  )J 

j     -g 2TT72 ~~^z  <9 

s         m    e 


where  the  Hamil ton-Jacobi  constant  coordinate,  which  normalizes  range  with 

respect  to  source  location,  is  zero  for  sources  located  on  the  z-axis.   It 

is  noted  that  the  integral  in  Eq.  (9)  is  a  line  integral  and  that  one  must 

integrate  along  the  branch  cut  and  around  any  encountered  branch  points 

(i.e.,  vertex  depths)  onto  the  next  Riemann  sheet.  The  node  in  Eq.  (3) 

at  z  =  *  for  trapped  rays  induces  C  ->  C  from  below  and  3C  /9C  ■*  1  from 

em  em 

below  as  z  ■*  ».     Thus  the  vertex  depth  for  the  ray  is  at  infinite  depth. 

Nevertheless,   the  range   quadrature,   Eq.    (9),   for  rigorous   ray  tracing 

remains  bound  and  integrable  for  trapped  rays  as  shown  in  Section   II. 

Furthermore   in  the  limit  C     ■*■  C     by  L'Hospital's  rule,   the  integrand  in 

em 

Eq.    (9)   goes   to  zero. 

For  cylindrical    symmetry  by  Eq.    (9),   the  ray  also   is  embedded  in  a 
radial-vertical    plane.     The  D/E  angle,   $,   for  the  ray  for  rigorous   ray 
tracing  by  Eq.    (9)    is   given  by 


4>    =     cot"1^) 

i      c   n  -  (C  /CJ3(3C  /3CJ] 
.-1  r     el        v  m'    e'         e       m     -.  /ln\ 

=     cot     {  p ?    -.  /? }  (10) 

[C  -  C  ]  u 
1  m    eJ 

which  in  general  differs  from  \p   as  given  by  Eq.  (8).  Hence  the  ray  direc- 
tion is  not  parallel  to  the  wave  normal.  There  is  a  precedent  for  rays  not 
being  normal  to  the  wave  surfaces  of  constant  phase  or  W  in  the  optics  of 


anisotropic   crystals.      In  a   preceeding  paragraph  of   this   section,    it   has 

been   shown   that   in   general    C     is  anisotropic. 

At   the   source,    3C  /3C     =   0   in   order  to   satisfy   the   boundary  condi- 
e       m 

tions   for   Eq.    (3).      Hence  at   the  source,   we   have   by  Eqs.    (8)   and    (10) 
that.  ip(z$)    =   $(zs)  . 

Equation    (10)    is   Snell's   law   for   the  ray   for   rigorous    ray  tracing. 
Because   Eq.    (3)    has  a   node  st  z  s  «,   we  have  that   v  ■*■  0  and  $  ■*  90°   as 
z  ■+  «>.     Hence  the  ray   is   parallel    to   the  wave  surface  at   the  vertex   point 
at  z   =  oo  for  rigorous   ray  tracing. 

It   is  noted  that  Ref.   1    needs  correction  because  Eq.    (8)  of  this  exposi- 
tion was   given  in  Ref.   1    as  the  applicable  modification  of  Snell's  law  without 
distinguishing   between  wave  normals  and  ray  directions.     Nevertheless  all 
applications   in  Ref.    1   of  Snell's   law,  which  were  based  upon  an  effective 
sound  velocity  profile  developed  only  to  first  order,  are  correct  because 

to  first  order  only,   C     is  a   function  of  both  z  and  co  but  not  C   .     Thus  a 

'       e  m 

first-order  C     is   isotropic.     We  note  that  if  C     is  independent  of  C   , 
e  em 

then  Eq.    (10)   reduces  to  Eq.    (8). 

We  offer  a  caution  with  regard  to  rigorous  ray  tracing  with  an  effec- 
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tive    sound  velocity  profile  as  given   by  Ce*(A-8)  where  $  is  the  solution 

to  Eq.    (3).     The  worker  is  admonished  to  use  C     only  in  a  Hamilton-Jacobi 
process   (regrettably,   it  does  not  render  correct  results  when  C     is  substi- 
tuted for  C  into  a  Lagrangian  or  Hamiltonian  calculus  of  variation  process). 


II.      ANALYSES   OF   SAMPLE   APPLICATIONS 
Let   us  assume  a   sound   velocity  profile,    C(z),    such   that 


C(2)      ■      [(1420F2  -  -0*0361   ,-V2i       2   >  0 

(1420) 


where  C  is  in  meters  per  second  and  z  is  in  meters.  The  source  depth  is 

assumed  to  be  zero.  This  profile  was  chosen  for  investigation  because 

(1)  it  represents  the  profile  that  is  dominated  by  the  pressure  gradient 

in  the  absence  of  a  thermocline  and  (2)  in  a  perturbation  expansion  the 

first  order  calculation  of  C  is  identical  to  C.  Here  we  examine  C  for 

e  e 

various  values  of  C     and  to.     Subsequently  we  rigorously  trace  a  particular 
ray  to  exhibit   its  penetration  of  the  classically  forbidden  region  beyond 
the  WKB  vertex  depth  and  to  exhibit  the  divergence  between  the  ray  direc- 
tion and  wave  normal. 

Let  us  now  numerically  solve  Eq.    (3)  and  exhibit  by  Fig.   1   the  effec- 
tive sound  velocity  profile,   C  ,  for  a  harmonic  surface  source  (z     =  0) 
with  to  =  50  rad/s   for  rays  with  initial    D/E  angles    [note  ip(z   )   =  <J»(z   ); 
z     =  0]  of  5°,  10°,   15°    (i.e.,   C     =  1425.4242  mps ,   1441.9058  mps,  1470.0922 
mps  respectively).     In  Fig.   1,   the  nodal    behavior  of  Eq.    (3)   is  manifested 

by  C     -*•  C     from  below  as  z  ■+  ».     Hence  we  observe  that  the  vertex  depth 
e        m 

recedes  to  z  =  ».     Well   above  the  WKB  vertex  depths,  C     and  C  are  good 
approximations  of  each  other  (this  may  not  be  true  for  less-well-behaved 
profiles   ),  and  by  the  initial   conditions     C     is  tangent  to  C  at  the  source 
depth  z    .     For  completeness,  the  actual   sound  velocity  profile,   C,   is  also 
exhibited  on  Fig.    1.     Note  that  in  the  limit    where   the  initial    D/E  angle 

goes  to  90°,   (i.e.,   C     ■*■«),  then  C     would  mimic  C  as  closely  as  possible 

1  /? 
for  the  value  of  A  (here  A  '     =  1420  m/s)   but  still   Cg  would  not  go  to  C. 


However  in  the  limit  A  ►  »  and  the  initial  D/E  angle  is  not  0°,  then  C   -  C 

(the  analogy  in  quantum  mechanics  would  be  the  infinite  energy  case). 

In  Fig.  2  we  exhibit  various  effective  sound  velocity  profiles  for 

surface  sources  with  an  initial  D/E  angle  of  10°  (i.e.,  C  =  1441  .9058  mps) 

and  with  w  =  5  rad/s,  50  rad/s,  500  rad/s.  The  sound  velocity  profile,  C, 

is  also  exhibited  on  Fig.  2.  We  note  that  for  finite  depths  above  the 

WKB  vertex  depth  C  approaches  C  as  w  increases  and  that  for  finite  depths 

below  the  WKB  vertex  depth,  C  approaches  C  as  u  increases.  In  the  limit 

r     e  m 

co  ■*  °°,  we  have  that 


c,   c  <  c 

—  m 
Lim  Ce  ->    f 

^  c  ,   c  >  c  . 

m        m 


Let  us  now  investigate  the  correspondence  principle  of  how  rigorous 
ray  tracing  goes  to  classical  ray  tracing  in  the  high  frequency  limit. 
In  the  case  co  -*■  °°  and  below  the  WKB  vertex  depth  (which  is  the  classical 
ray  tracing  vertex  depth)  in  the  classical  forbidden  region,  the  ray  direc- 
tion is  parallel  to  the  wave  front  and  so  the  ray  is  embedded  in  a  wave 
surface  of  constant  phase.  Hence  below  the  WKB  vertex  point  for  co  ■*■  », 
the  ray  senses  a  constant  phase  and  thus  the  ansatz,  Eq.  (2),  in  the  forbidden 
region  has  a  constant  phase  despite  to  -►  °°,  that  is,  the  integrand  of  Eq. 
(4)  behaves  such  that  (A  -  C   -  $)  -*■  0  as  co  ■*  «.  The  boundary  condition 
for  infinite  frequency  for  the  solution  to  Eq.  (1)  requires  that  in  the 
classical  forbidden  region  that  e(z)  =  0  for  z  below  the  WKB  vertex  depth. 
Therefore  we  have  that 


(6+  -  6_)/2, 


CO 


10 


(* 


lOl      Wipthd     :      low    the     ..I..'     '/.■:•      x    -.,-'.,tl 


1 

which  goes   to   zero  by  L 'Hospital ' s   rule  despite  the  singularity   in   the 

amplitude    factor   in   both   0+  and  9     as   given   by   Eq.    (2).      Thus   for  u  ■+  «>, 

there   is   no   penetration  below  the  WKB   vertex   depth,   and  by  Eq.    (11)    below 

the  WKB  vertex   depth  the   integrand  of  the  range  quadrature,   Eq .    (9),   for 

rigoreus   ray  tracing   is   zero   by  L'Hospital's   rule  while  above  the  WKB 

vertex   point   Eq.    (9)   reduces  to  the  range  quadrature,   Eq.    (A5),   for  classical    ray 

tracing  as   3C  /3C     ->  0.     Thus   in   the   limit  w  -*•  °°,   rigorous   ray  tracing   corres- 

3 
ponds  to  classical    ray  tracing   [the  quantum  mechanical    analogy     is   the  corres- 

o 
pondence  principle  for   (h  /2u)  ■*■  0  where  R  is  Planck's  constant  and  y  is  mass]. 

We  exhibit  on   Fig.   3  the  rigorous  ray  trace  for  a   surface  source  for 
co  =  50  rad/s  and  the  <j>(z    )   =  ip(z   )    =  10°    (i.e.,   C     =  1441.9058  mps)   by 
numerically  evaluating  Eq.   (9).     The  corresponding  classical   ray  trace  is 
exhibited  also  on   Fig.   3  for  a   half-cycle  in  range  (i.e.,  until    the  first 
WKB  vertex  depth  is  encountered  at  z  =  1189.1117  m)   where  classically 
r  =  13,490.938  m.     While  the  corresponding  rigorous  ray  trace  has   its 
vertex  depth  at  z  =  °°  which  is  not  exhibited,  we  report  that  the  range  for 
the  corresponding  rigorous   ray  trace  to  the  first  vertex  depth  has  been 
evaluated  to  be  13,557  ±  6  m.     This  infers  that  classical   ray  tracing  for 
one-half  cycle  for  these  conditions  has  an  estimated  error  of  65  m  in  range 
or  about  0.498%  in  range. 

We  note  that  while  the  vertex  depth  for  rigorous  ray  tracing  has  re- 
ceded to  infinity,   the  range  quadrature  has  remained  finite.     We  also  report 
that  the  time  quadrature,  which  is  the  action  variable  for  a  full   cycle 
remains   finite  as  inferred  by  the  analogous  quantum  mechanical   effective 

action  variable  remaining  finite  for  bound  states  of  the  symmetric  bi-linear 

3  4 
potential   whose  eigenfunctions  are  offset  Airy  functions.' 

On   Fig.   4  we  exhibit  as  a  function  of  depth  the  wave  normal,  \\jt  and 

the  ray  angle,   <j>,   for  the  ray  that  was  rigorously  traced  on  Fig.   3.     For 
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classical  ray   tracing,  which  also  is  exhibited  on  Fig.  4,  the  wave  normal 
and  ray  angle  are  equivalent.  We  note  that  well  above  the  WKB  vertex  depth 
in  the  classically  allowed  region  both  \p   and  0  approach  their  classical  ray- 
tracing  prediction.  But  below  the  WKB  vertex  depth  we  have  that  as  2  ■*•  », 
then  ip  ->  0  and  4>  ■*  90°.  Hence,  as  z  ■>»,  the  ray  travels  parallel  to  the 
wave  front  as  noted  earlier. 
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APPENDIX  A 
CLASSICAL   RAY   TRACING 

In   this  appendix,   we  outline   the  development  of   the  Hamil ton-Jacobi 
theory  of  classical    ray  tracing   in  the  short  wavelength,   A,  or  high  frequency 
1  imit   (i.e.,  A  ->  0  or  co  ■*  °°)   and  show  that  the  ray  path   is  always  normal    to 
the  wave  fronts.     Here  we  assume  that  the  sound  velocity  profile,   C,   is 
only  dependent  upon  depth,   z,  that  the  earth  is   flat  and  that  the  ocean 
is   bottomless.     We  use  cylindrical    coordinates  and  assume  the  source  is 
located  on  the  z-axis  at  depth  z    . 

We  start  with  the  Lagrangian,   L,   for  ray  tracing  that  for  azimuthal 
symmetry  is  given   in  cyclindrical    coordinates   by 


L(zz)     -     H+i2)1/2 


where  z  =  dz/dr  where  r  is  the  radial  coordinate  and  for  sources  on  the  z- 
axis  r  determines  range.  The  conjugate  momentum,  P  ,  with  respect  to 
depth  is  given  by 


p,    =   -Ml)  sin 


z        »z    '    '    c(z)(i+i2),/2  «* 


where  <f>  is  the  depression-elevation  (D/E)  angle  of  the  ray  (for  azimuthal 
symmetry,  the  ray  is  embedded  in  the  radial-vertical  plane).  As  L  is  not 
an  explicit   function  of  range,  r,  the  classical   ray-tracing  Hamiltonian,  H, 
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is  given  by 

H(z,Pz)  =  zPz  -  L  =  -[C2   -  P2z)]/2 

cos  <$>  1 

"  CTz)      "  C 
m 


(AT) 


which  gives  Snell's  law  since  H  is  a  constant  of  the  motion  for  H  indepen- 
dent of  r.  The  constant  of  the  motion  for  any  particular  ray  may  be 
expressed  in  terms  of  its  vertex  velocity,  C  .  The  Hamil ton-Jacobi  equa- 
tion for  classical  ray  tracing  is  given  by 


3S/3r  +  H(z,3S/3z)   =  0 


or 


3S/3r  -    [C(z)"2   -   (3S/3z)2]1/2     =     0  (A2) 

where  S   is  Hamilton's   principal    function  for  classical   ray  tracing  and  has 
the  dimensionality  of  time.     Thus  the  surface  of  constant  S  represents  a 
wave  surface.     Because  the  only  range  dependence  is  restricted  to  the 
3S/3r  term   in  Eq.   (A2),  we  may  assume  a   solution  of  the  form 


S(z,Cm,r)      =     W(z,Cm)   +  r/Cm  (A3) 


where  W   is  Hamilton's  characteristic   function  for  classical    ray  tracing 
From  Eqs.    (A2)  and   (A3),  W  may  be  expressed  as 
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W      =     /      [C~2   -   C~2]1/2   dz.  (A4) 


The  wave  normal    is   given  by  the  direction  of  the   gradient  of  Hamilton's 
principal    function,   S,  which  by  Eqs .    (Al)   and   (A2)    is   given   for  azimuthal 
symmetry  as 

1     C  1  11/2" 

vs  =  r  V  +  [J2  -  V     V 

m     r         C         C 
m 

For  azimuthal    symmetry  the  wave  normal    is   embedded   in  a   radial-vertical 
plane  and   the  wave-normal   may  be  described  by  the   D/E  angle,  ^,   of  the 
wave  normal    as  given  by 

C2 
tan"1^-  lf1/2     =    4, 

r 

But  the  equation  of  motion  for  the  ray  for  classical    ray  tracing  from  a 
source  on  the  z-axis   is  given  for  azimuthal    symmetry  by 


9S                    n                  (Z        'c   dz'  t   r  /*K\ 

-     0     -        j     — p PH72  +  r  (A5) 


^W  z;      [C^C2}1'2 

s  m  J 


where  the  z   integral    is  a  line  integral   with  branch  points.     For  azimuthal 
symmetry,   the  ray  is  also  embedded  in  the  radial-vertical    plane  and  the 
ray's  D/E  angle,  <j>,   is  given  by 

C2 
<j>     =     tan'^dr/dz)"1      =     tan[^-l]1/2     =     ^. 
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Hence  for  classical  ray  tracing,  the  D/E  angle  for  both  the  ray   and  the 
wave  normal  are  equal  to  each  other.  Therefore  by  azimuthal  synmetry, 
the  ray  path  is  always  normal  to  the  wave  fronts. 
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FIGURE    CAPTIONS 

1.  The  effective   sound  velocity  profiles    for  a   harmonic   source  at  zero 
depth  with  uj  =  50  rad/s   for   initial    D/E  angles  of  5°,   10°  and  15° 
(i.e.,   Cm  =   1425.4242  mps,   1441.9058  mps   and  1470.0922  mps   respec- 
tively).    The  sound  velocity  profile,    C,   is   represented  by  the 
dashed  curve. 

2.  The  effective  sound  velocity  profiles  for  a  harmonic  source  at  zero 
depth  with  an  initial    D/E  angle  of  10°   (i.e.,   C     =  1441 .9058  mps) 
for  oo  =  5  rad/s,   50  rad/s  and  500  rad/s.     The  sound  velocity  profile, 
C,   is   represented  by  the  dashed  curve. 

3.  The  rigorous  ray  trace  for  a  harmonic  sound  source  at  zero  depth,  with 
oj  =  50  rad/sec,  with  an  initial    D/E  angle  of  10°  for  the  sound 
velocity,   C,  exhibited  on  Figs.   1   and  2.     The  classical   ray  trace  is 
represented  by  the  dashed  curve. 

4.  The  deflection/elevation  angle  for  the  ray,  <j>,  and  wave  normal,  ty, 
as  a  function  of  depth  for  the  rigorous  ray  trace  exhibited  by 
Fig.   3.     The  deflection/elevation  angle  for  classical   ray  trace  is 
exhibited  by  the  dashed  curve. 
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